Yukawa Institute Kyoto 



DPSU-11-2 
YITP-11-18 



The Exceptional (Xi) (Q')-Racah Polynomials 



o 

(N 



> 



Satoru Odake" and Ryu Sasaki^ 



Department of Physics, Shinshu University, 



■ Matsumoto 390-8621, Japan 



Yukawa Institute for Theoretical Physics, 
Kyoto University, Kyoto 606-8502, Japan 



Abstract 



The exceptional Racah and g-Racah polynomials are constructed. Together with 
the exceptional Laguerre, Jacobi, Wilson and Askey- Wilson polynomials discovered by 
fS| ' the present authors in 2009, they exhaust the generic exceptional orthogonal polyno- 

mials of a single variable. 
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§ : 1 Introduction 

The exceptional (X^) (g)-Racah polynomials and related exceptional orthogonal polynomials 
' constructed as the main part of the eigenf unctions of the shape invariant and exactly 

[ solvable discrete quantum mechanics with real shifts [1], which are deformations of those 
governing the corresponding orthogonal polynomials, i.e. the (g)-Racah polynomials, etc. 
[21 [31111 [5]. The method of deformations is essentially the same as that for the (Xi) Wilson and 
Askey- Wilson polynomials derived by the present authors in 2009 [6j. Namely, the potential 
functions of the original Hamiltonians are multiplicatively deformed in terms of a degree i 
eigenpolynomial with twisted parameters. The exceptional (g)-Racah polynomials and the 
exceptional Wilson and Askey- Wilson polynomials share many properties. One pronounced 
difference is that there are only finitely many exceptional (g)-Racah polynomials in contrast 
with the infinitely many types of the exceptional Wilson and Askey- Wilson polynomials. 
For example, starting from the (g)-Racah polynomials of the highest degree N, there exist 
N — 1 different types of the exceptional (g)-Racah polynomials, for which the highest degree 



is always N. On the other hand, there are infinitely many different types of the exceptional 
little g-Jacobi polynomials, since the degrees of the original little g-Jacobi polynomials are 
not bounded. These exceptional (X^) polynomials are exceptional in the sense that they 
form a complete set of orthogonal polynomials in spite of the fact that the lowest member 
of the polynomials has degree (> 1) instead of a constant. Thus they do not satisfy the 
three term recurrence relations. 

Historically the Xi Laguerre and Jacobi polynomials were discovered by Gomez-Ullate et 
al [7] in 2008 within the framework of the Sturm-Liouville theory. Soon they were rederived 
as the main part of the eigenf unctions of shape invariant quantum mechanical Hamiltonians 
by Quesne and collaborators |8j. In 2009 the present authors derived the infinitely many 
Xi Laguerre and Jacobi polynomials by deforming the Hamiltonian systems of the radial 
oscillator and the Poschl- Teller potential in terms of the eigenpolynomials of degree i [9l 
[Tot [TT] . The examples of Gomez-Ullate et al and Quesne et al are the first members of 
the infinitely many exceptional polynomials. For the recent developments of the exceptional 
orthogonal polynomials, see fl2\ \T3\ [HI [151 [IS]- K is worth remarking that the general 
knowledge of the solution spaces of exactly solvable (discrete) quantum mechanical systems 
governed by Crum's theorem [17] and its modifications [HI [191 [20l [21] has been very helpful 
for the discovery of various exceptional orthogonal polynomials. 

The orthogonal polynomials of a discrete variable [2] have played important roles in 
many disciplines of physics and mathematics [2l |3l [1] . See [22] for recent applications. Let 
us comment on the birth and death processes, the typical examples of Markov chains, which 
could be considered as a discrete version of the Fokker-Planck equations [23j. As shown in 
[21], the explicit examples of 18 orthogonal polynomials in [Tj, the (g)-Racah, (g)-(dual)- 
Hahn etc, provide exactly solvable birth and death processes [H [25] . That is, for the given 
birth and death rates {B{x), D{x)} which define the Hamiltonian (12. ip . the corresponding 
transition probabilities are given explicitly, not in a general spectral representation form of 
Karlin-McGregor [26j . The exceptional versions presented here also provides ample examples 
of exactly solvable birth and death processes. 

The present paper is organised as follows. In section two, the basic principles of the 
discrete quantum mechanics with real shifts are briefiy reviewed with an emphasis on the 
shape invariance. The details of the Racah and g-Racah polynomials are recapitulated in 
section three. The exceptional Racah and g-Racah polynomials are introduced in section 
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four. The intertwining relations connecting the original (g)-Racah and the exceptional (g)- 
Racah polynomials are explored in section five. These two sections are the main part of this 
paper. Several exceptional orthogonal polynomials, the dual (g)-Hahn and the little g-Jacobi 
polynomials are derived from the exceptional (g)-Racah polynomials in section six through 
certain limiting processes. The final section is for a summary and comments. 

2 General Setting: shape invariance 

Let us recapitulate the essence of the discrete quantum mechanics with real shifts developed 
in [l]. The Hamiltonian "H = {'Hx,y) is a tridiagonal real symmetric (Jacobi) matrix and its 
rows and columns are indexed by non-negative integers x and y, x,y = 0,1, . . . , x^ax, either 
finite (xmax = N) or infinite (xmax = oo). The Hamiltonian H has a form 

n = ~^/^^e^^/D{^- ^/I^e-^ (2.1) 
n.^y = -^/B{x)D{x + l)6,+i,y - y/B{x-l)D{x)6,^i,y + {B{x) + D{x))5x,y, (2.2) 

in which the two functions B{x) and D{x) are real and positive but vanish at the boundary: 

B{x) > 0, D{x) > 0, D{0) = ; B{xme.x) = for finite case. (2.3) 

The Schrodinger equation is the eigenvalue problem for a hermitian matrix (^max = N or 
oo), 

'H4>n{x) = £n4>n{x) (n = 0, 1, . . . , nmax), = £q < £i < £2 < ' ' ' ■ (2.4) 
The Hamiltonian (12.11) can be expressed in a factorised form: 

T~i jA., {^^x,y)-i ^'^y,x)i (-^j V 0, 1, ... , X^iax); (2"^) 

^ = /BO^-e^/D(^), = ^B{x)-^D{x)e-^, (2.6) 
Ax,y = \/B{x) 5x,y - V-D(x + 1) 6x+i,y, {A^)x,y = \/B{x) 5x,y - V-D(x) S^-i^y (2.7) 

The zero mode A(t)o{x) = {(poi^x) > 0) is easily obtained: 0o(a^)^ = 11^=0 i^y+i) ' "^i^h 
the normalization 0o(O) = 1 (convention: nfc=n * ~ -'-)• adopt the standard euclidean 
inner product ( , ) of two real functions on the grid as {f,g) ^!^=o' /(^)fl'(^)- Then the 
orthogonality relation reads 

{4>n, (Pm) = Snm (n, m = 0, 1, . . . , n^ax) • (2.8) 
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Here is the normalization constant to be specified later. 

Shape invariance, a sufficient condition for the exact solvability [H [23 [281 129] , is realised 
by specific dependence of the potential functions on a set of parameters A = (Ai, A2, . . .), 
to be denoted by B{x; A), D{x; A), ^(A), 'H(A), £^„(A), 0„(a;; A) etc. The shape invariance 
condition is 

A{X)A{Xy = kA{X + 5)U(A + 6)+ £i{X), (2.9) 

where 5 is a certain shift of parameters and /t is a positive constant. It should be stressed 
that the above definition is much stronger than the original definition by Gendenshtein |30j . 
The shape invariance condition (12. 9p combined with the Crum's theorem pTl [T9l l20l [2T| 
implies that the entire energy spectrum and the excited states eigenfunctions are expressed 
in terms of Si{X) and (poix; A) as follows: 

71-1 

^n(A) = 5^fi:%(A + s(5), (2.10) 

s=0 

(t)n{x] X) oc ^(A)U(A + (5)U(A + 25)^ ■■■A{X + {n- 1)5)^ ^x] X + n5). (2.11) 



^(A)0„(x; A) = -j=== fn{X)<Pn~i {x; X + 5), (2.12) 
±5(0; A) 



We have also 



^(A)Vn-i(x;A + 5) = /b(o;A)6„_i(A)0„(x;A), (2.13) 

where /n(A) and 6„_i(A) are the factors of the energy eigenvalue, £n{X) = /„(A)6„_i(A). 

For the (g)-Racah and the other polynomials to be discussed in the present paper, the 
eigenfunction has the following factorised form, 

0„(x; A) = ^)Pn{x- A), P„(a:; A) = P^x', A); A), (2.14) 

where Pn{i]{x; A); A) is a polynomial of degree n in the sinusoidal coordinate ri{x; A). The 
sinusoidal coordinate considered here is a monotone increasing function of x satisfying the 
boundary condition ?7(0; A) = [H [29]. We choose the normalization 

P„(0;A) = 1, (2.15) 

and set P_i(x; A) = 0. For later convenience, let us remark on the relation 

PMl; A); A) = P„(l; A) = 1 - (2.16) 



The orthogonality relation (12 .Sp becomes 

^ 0o(x; XfPnjx; X)Pmix; A) = 5nm (n,m = 0, 1, . . . ,nmax)- (2.17) 



x=0 



The forward shift operator J-'(A) = (J^J;^y(A)), the backward shift operator B{X) = {Bx,y{X)) 
and the similarity transformed Hamiltonian 7^(A) = {T-Lx.yW) i^yV — 0, 1, . . . , Xmax) are 
defined by 

J-(A) = V5(0;A) 0o(a;; A + 5)"^ o ^(A) o 0o(x; A) 

= 5(0;A)y.(x;A)-i(l-e^), (2.18) 
i3(A) 0o(x; A)-i o ^(A)t o 0o(x; A + 5) 

V -°(^; 

^ (5(x; A) -D(x; A)e"^)¥?(x; A), (2.19) 



5(0; A 

•H(A) =^ 0o(a;; A)-i o •H(A) o 0o(a;; A) = i3(A)J^(A) 

= B{x; A)(l - e^) + Z?(x; A)(l - e"^), (2.20) 



where the auxiliary functions ip{x) is defined by pQ 



def B{0;X) (f)o{x;X + 6) 7]{x + 1; X) - 7]{x; X) 

^{x, X) = J — — TT— -7 — = 7--rT , ^(0, A) = 1. (2.21) 

\l B{x;X} 0o(a;;A) ^7(1; A) 

Their action on the polynomials is {n = 0,1, . . . , rimax) 

^(A)P„(x; A) = /„(A)P„_i(x; X + S), (2.22) 

S(A)P,_i(x; A + (5) = 6„_i(A)P„(x; A), (2.23) 

n{X)Pn{x; A) = ^„(A)P„(x; A). (2.24) 

The above difference equation fl2.24p for the polynomial P„ reads explicitly as 

B{x) (Pnivix)) - Pnivix + 1))) + D{x) (P„(r/(x)) - P„Mx - 1))) = SnPnMx)), (2.25) 

in which the parameter dependence is suppressed for simplicity. 

3 Original Systems: (Q')-Racah polynomials 

Here we present various properties of the Racah (R) and the g-Racah (gR) polynomials as 
explored in [1]. In general there are four cases of possible parameter choices indexed by 
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(e, e') = (±1,±1). Here we restrict ourselves to the (e, e') = (1,1) case for simplicity of 
presentation. 

The set of parameters A, which is different from the standard one (a, /3,7, 6) [5j, its shift 
5 and K are 

R: X = {a,b,c,d), (5 = (1, 1, 1, 1), k = I, (3.1) 
qR: = {a,b,c,d), (5 = (1, 1, 1, 1), k = q~\ < q < 1, (3.2) 

where q^ stands for g(-^i'-^2' - ) = (g-^i, q^'^^ . . .). We introduce a new parameter d defined by 

7defJa + 6 + c- (i-l:R , , 

~ \ ahcd-\-^ : qR ' ^^'^^ 

The Hamiltonian is a finite dimensional matrix and the maximal values of x and n are 
a^max = ''^max = N and we could choose 



R : a = -A^ or h = -N or c = -iV, 

gR : a = or b = q^^ or c = (3.4) 



to ensures the boundary condition for B, B{xma.x) = 0. The potential functions B{x] A) and 
D{x; A) are 

{x + a){x + b) {x + c){x + d) _ ^ 

R( \\ - ) {2x + d){2x + 1 + d) ■ , . 

: qR 



(l-dg2^)(l-rfg2x+i) 
{x + d — a){x + d — b){x + d — c)x 



R 



n( w - ) {2x - 1 + d){2x + d) ' , . 

~ \ ~{l-a-^dq^){l-b-^dq'=){l-c~^dq^){l-q'') ' ^^'''^ 

[ (l-rfg2-i)(l_rfg2x) • 

The parameter ranges are restricted by the positivity of B{x; A) and D{x; A). When we need 
to specify them, we adopt the following choice of the parameter ranges: 

R : a = -A^, a + b> d>0, < c < 1 + d, 

qR : a = g"^, < ab < d < 1, qd < c < 1. (3.7) 

The energy eigenvalue and the sinusoidal coordinate are 

p /\\ _ j T^i^^ + d) :R , ,._(x{x + d) :R , . 

^n[^) - I ^g-n J^„) . ' V[x, Aj - I _ ^^^^ _ _ . (3.8j 
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The eigenfunctions have the factorised form fl2.14p and the orthogonal polynomials are the 
Racah and the g-Racah polynomials: 



P„(x; A) =P„(r/(x;A);A) 



4-^3 



— n, n + d, —X, X + d 



a, b, c 
dg", q~^, dq 

a, b, c 

Rn{i]{x; A); a — 1, d — a, c — 1, d — c 



4<P3 



Ij : R 

q;q) : qR 



R 



(3.9) 



(3.10) 



\ Rn{l + d + r](x; X); aq ^,da ^,cq ^,dc ^\q) : qR 
Here -R„(- ■ ■ ) are the standard notation in [S]. The auxiliary function <f{x; A) (I2.2ip reads 

2x + d+l 

ip{x; A) 



d + 1 



■ R 



(3.11) 



1 — dq 



: gR 



The constants /n(A) and 6ri(A) appearing in f l2.12l) - fl2.13l) are 

/„(A)=^„(A), 6„(A) = 1. 

The orthogonality measure (polx; A)^ and the normalisation constants (i„(A)^ are 

(a, b, c, d)x 2x + d 

I 7i ^ 1 

(x; A)2 



(3.12) 



cin(A)^ 



(1 + c/ — a, l + c? — 6, 1 + c? — c, 1)^ d 
(a, b,c,d] q)x 1 — (ig^^ 

{a~^dq,b~^dq,c~^dq,q;q)xd^ 1 — d 

{a,b,c,d)n 2n + d 

{1 + d — a, 1 + d — b, 1 + d — c, l)n d 

1 + d — a,l + d — b,l + d — c)n 



R 
qR 



(3.13) 



X ■ 



{d+l)N{d+l)2N 

[a, b,c,d; q)n 1 — dq^"^ 



: R 



(3.14) 



(a ^dg, b~^dq, c~^dq, q ; g)„ 1 — 



\7V/ 



-irfg, c-^c/g ; g),v rf^g^^(^+^) 



{dq;q)N{dq;q) 



: gR 



2Ar 



4 Deformed Systems: (Q')-Racah polynomials 

For each positive integer £=1,2,...,A^ — 1, we can construct a shape invariant system by 
deforming the original system {£ = 0) in terms of a degree i eigenpolynomial ^£ of twisted 
parameters. 
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We set 



^L. = N-i, nL^'^'N-i, (4.1) 



and take 

R : a = -A^ or b= -N, 

qR: a = q-^ or b = q'^. (4.2) 

The deforming polynomial which is a polynomial of degree £ in r]{x; A + (£ — 1)6), is 
defined from the eigenpolynomial Pe{x; A): 

ieix;X)=Uvix;X+ii-l)S);X) 



def 



Pe{x;i{X + {£-l)S)) : R,gR 

■£, i — a — b + c + d — 1, —X, x + d + 



(i — a, (i — 6, c + ^— 1 



'^'^ a-^d, b-^d, cq'-^ 



Ij : R 
g ; g ) : gR 



(4.3) 



which satisfies the normalization 

6(0; A) = 1. (4.4) 
Here the twist operator t acting on the set of parameters A = (Ai, A2, A3, A4) is 

t(A) = (A4 - Ai, A4 - A2, A3, A4) : R, gR. (4.5) 

This is the most important ingredient of the deformation. For the appropriate parameter 
ranges, for example as given in 03.71) . the deforming polynomial C,e{x] X) is positive at integer 
points X = 0, 1, . . . ,x^^^ + 1, because the polynomial ^i{y; X) has no zeros in the interval 
< y < vi^max + I; X + {£ — 1)6). It satisfies the following two formulas, which will play 
important roles in the derivation of various results: 
1 



ip{x; X + iS + S) 



1 



(v^ix; X + i6)- v^ix; X + M)e^)6(x; A) = Ao(A)6(a;; A + S), (4.6) 

. X ^^-^U{x;X+{i-l)S)-v^{x;X+{i-l)6)e-^)Ux;X + 6) 

Lp{x; X+ {i - 1)0 + 0) ^ ^ 

= be,o{X)Ux;X). (4.7) 

Here ff (x; A), (x; A), ff (x; A), f|'(x; A) are the factors of the potential functions B{x; X) 
and D{x; A): 

d-^{x + a){x + b) :R 
<(.;A,-^ (4.8) 

1 — d 



where 6 is 



d~^{x + c){x + d) :R 
vn..X) = { . (4.9) 

1 — d 

d~^{x + d — a)(x + d — b) : R 

(/"-"^(a; + (i — c)x : R 

.-(.;A)-^ (4.11) 

B<r,X)--^. . (4.12) 
A + d)(y9(x + i; A + d) 

D{x; A) = -v^ ^ ^ J ^ ^ \' ' ^, (4.13) 

A + d)(f{x - j;X + d) 



5 = (0,0,-1,-1) : R, qR. (4.14) 



The constants fi^nW and be^nW are given by 

, , , ,c + 2i + n-l ^ 
I [a + b-d + n) — ; :R r , i o 

g-(l-a6ci-V) \_^^^,_, :gR ^ ^ " ' 



The eqs. fl4.6p - fl4.7p are identities relating ^e{x] A) and C,e{x; X + S). They are reduced to the 
identities satisfied by the (basic) hypergeometric functions, (2.74)-(2.75) in [15]. Note that 
these two equations (14. 6p - (14. 71) imply the difference equation for the deforming polynomial, 

t(A + 1)6)) (1 - e^) + D{x; t(A + (£ - 1)S)){1 - e-^))6(x; A) = £^{t{X))Ux; A), 

(4.16) 

which corresponds to (I2.24p . 

Let us introduce new potential functions Bi{x; X) and Dg^x; X) by multiplicatively de- 
forming the original ones in terms of the polynomial ^^(x; A): 

B,ix; A) B{x- X + i5) J'^^'^^^ + I; X + S) 

D,{x- A) D{x- A + ^s^ Ux+^-^>^)Ux-l;X + S) _ 

it{x] X) ii{x] X + 5) 

See the corresponding expressions for the exceptional Wilson and Askey- Wilson polynomials 
(30)-(31) of [6j and (2.42)-(2.43) of [15]. They define a deformed Hamiltonian U^iX) = 
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{'He;x,yW) and other operators A (A) = {Ai-x,yW) and A(A)^ = {{Ae{Xy)x,y) = (^<?;y,x(A)) 
{x,y = 0,l, ...,a;^ax) by 

= A(A)U,(A), (4.19) 
MX) = y/Be{x; A) - eVA(x; A), ^,(A)^ = ^^.(x; A) - ^/D,{x;X) e'^ . (4.20) 

We have -0^(0; A) = and -^^(x^ax) ^) — 0- The parameter ranges are restricted by the 
positivity of B^{x] A) and Di{x] A). When we need to specify them, we consider the parameter 
ranges fl3.7p . 

The deformed system is shape invariant, too: 

A(A) A(A)^ = kA,{\ + 5)U,(A + 5)+ £,,i(A), (4.21) 

or equivalently, 

^Bi{x+l]\)De{x + l-X) = K^/Be{x; A + S)Di{x + 1; A + 6), (4.22) 

Bi{x; A) + Di{x + 1; A) = k{Bi{x; X + S)+ Di{x; X + S))+ ^£,i(A). (4.23) 

The proof is straightforward by direct calculation. In order to verify (14.231) . use is made of 
the two properties of the deforming polynomial C,e{x] A) fl4.6l) - fl4.7l) . 

The Schrodinger equation of the modified system is (n = 0, 1, . . . , n^ax) 

H,(A)0,,„(x; A) =^,,„(A)0,,„(x;A), ^,,„(A) = f„(A + M). (4.24) 

The ground state 4>efl{x; A), which is annihilated by ^^(A), is 



>; A) 



\fiDe{y + l;X) 



Mx; A) 0o(x; A + iS) J ^'f'l^l, (4.26) 

V X)Oix + 1; X) 

with the normalisation 0^,o(O; A) = 1 and ^/'^(O; A) = 1. The excited states wavefunctions 
have the factorised form as fl2.14p : 

0^,„(x; A) = M^-^ X)Pi,n{x; A). (4.27) 

The exceptional (X^) {q)-Racah polynomial Pe^n{x', A) is bilinear in the deforming polynomial 
and the original polynomial P„: 

Pe,n{x;X) = PMx;X + iSy,X) 
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def 



1 



1 




This is one of the main results of the present paper to be compared with the similar expres- 
sions for the exceptional Laguerre & Jacobi polynomials (2.1)-(2.4) in [12], (2. 31), (2. 33) & 
(3. 37), (3. 40) of [H], for the exceptional Wilson & Askey- Wilson polynomials (2.52) in [T5] . 
The overall multiplicative factor is so chosen and as to realise the normalisation condition 



which can be shown by using (12.161) . This is a polynomial of degree i + n in ri{x; X + iS) . 
Note that Pifi{x; X) = ^i{x; X + 6) due to (14. 6p . which is obviously a polynomial of degree 
i in ri{x; X + iS). The exceptional orthogonal polynomial Pe^niu] X) has n real zeros in the 
interval < y < ri{x^^^; X + £d) for the appropriate parameter ranges, for example the range 
(13. 7p . It has i extra zeros which are usually complex and lie outside the above interval. 
The action of the operators Ai{X) and Ai{Xy on the eigenfunctions is 



Like the corresponding formulas of the original systems (I2.12p -( l2.13p . these are simple con- 
sequences of the shape invariance and the normalisation of the eigenfunctions. In the next 
section, we will derive these formulas through the intertwining relations and without re- 
course to the shape invariance of the deformed system (14.211) . The forward shift operator 
J-'e{X) = {J'i-x,y{X)), the backward shift operator ^^(A) = {Be-x,yiX)) and the similarity 
transformed Hamiltonian Tie^X) = {T-it,x,y{,X)) (x, ?/ = 0, 1, . . . , x^ax) are defined by 



P,,„(0; A) = 1, 



(4.29) 




(4.30) 



Ai{X)^i^n-i{x]X + 8) = a/ Bi{0] A) 6£,„_i(A)</),,„(a;; A), 
fe,nW = fn{X + id), 6,,„-i(A) = 6„_i(A + £6). 



(4.31) 
(4.32) 



J-,(A) = ^/B^iO^M^; X + 6)'' o AiiX) o i;,{x; A) 




— - fe(x + 1; A + (5) - ie{x; X + d)e' 
1, Aj V 



) 



(4.33) 




1 1 



B{0;X + i5)^^{x;X + S) 
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X (^B{x; X + i6)^i{x; A) - D{x] X + i6)^i{x + 1; X)e~^^<^{x; X + £5), (4.34) 
-HeiX) A)-' o -HeiX) o ij,{x- X) = B,{X)J',{X) 

= B{x;X + iS) Jl^^l^ rUx + i;X + S) _ 
'Ux + l-X)\ U^;X + 6) 

Xfix + l]X) fUx-l-X + 5) _Q 



+ D{x- X + i6 r\ ' ' . - e-j . (4.35) 
^e{x;X) V ^f,{x;X + d) J 

Compare with the similar expressions for the Laguerre & Jacobi polynomials (3.2)-(3.5) 

in [12], and for the Xji Wilson & Askey- Wilson polynomials (2.58)-(2.63) in [15]. Their 
action on the polynomials is (n = 0, 1, . . . , n^ax) 

J-,(A)P,,„(x; A) = /,,„(A)P,,„_i(a;; A + 5), (4.36) 

Bt{X)Pt,n-x{x- A + 5) = 6,,„_i(A)P,,„(x; A), (4.37) 

'Ht{X)P,,n{x- X) = £i,n{X)Pi,n{x; A), £e,n{X) = £niX + id). (4.38) 

The orthogonality relation is 

^tf^P,^4x;X)PU^-,X) = -^^ {n,m = 0,1,..., J. (4.39) 
The normalisation constants de^n{Xy are 

be,nW Si{X) be,n{X) /o,„(A + £6) s,{X) 

where S£{X) is defined by 

{d — a){d — h) 



: R 



^^(^) = , (i-a-M)(i-rM) ^ • (^-^^^ 

This will be proved in the next section. In the second equality of f l4.40p use is made of the 
explicit forms of (i„(A)^ fl3.14p . Note the positivity of the quantities, /£^„(A), 6^,„(A), ^^(A) > 
0. 

5 Intertwining Relations 

Here we demonstrate that the Hamiltonian systems of the original polynomials reviewed 
in §|3] and the deformation summarised in §|1] are intertwined by a discrete version of the 
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Darboux-Crum transformation. This provides simple expressions of the eigenfunctions of the 
deformed systems fl4.28p in terms of those of the original system, which is exactly solvable. 
It also delivers a simple proof of the shape invariance of the deformed system. The line of 
arguments goes parallel with those for the other exceptional orthogonal polynomials [El [15] . 

First let us discuss the general scheme. For an adjoint pair of well-defined operators 
^^(A) and ^^(A)^, let us define a pair of Hamiltonians ^^^^(A) 

Hf\\) = i,(A)U,(A), V!--\X) ''^ Ae{\)Ae{X)\ (5.1) 

and consider their Schrodinger equations, that is, the eigenvalue problems: 

?if)(A)0fj(x;A) =4t^(A)0fJ(x;A) (n = 0, 1, 2, . . .). (5.2) 

Obviously the pair of Hamiltonians are intertwined: 

?i^i+\A)i,(A)^ = A,{\)^Ae{X)A,{Xy = A,{\)%-\X), (5.3) 
A,{X)n^-^\\) = ^(A)i,(A)U,(A) = H\-\x)Ai{\). (5.4) 

If Ai{\)(j))^^{x]\) 7^ and Ai{\y (jy'f^ [x] X) ^ 0, then the two systems are exactly iso- 
spectral and there is one-to-one correspondence between the eigenfunctions: 

4^^W=4n(A), (5.5) 
05;) (x; A) oc ^(A)0(+J(x; A), 0{+J(x; A) (x A,{\)^ ^^^[^{x- A). (5.6) 

In the following we will present the explicit forms of the operators ^^(A) and ^^^(A)^, 
which intertwine the original systems in §[3] and the deformed systems in §|4l The operators 
At{\) = (A;x,2y(A)) and Ai{Xy = i{Ae{Xy)^,y) = iAe;y,xW) i^^V = 0, 1, • • • , x^ax) are 
defined by 



^(A) '=^' ^jBeix; A) - A(x; A), i,(A)t = ^J B^ix; A) - ^J D^ix; A) e'^ (5.7) 
where Bi{x; A) and De{x] A) are given by 

Mx; A) B (x; t(A 1)5)) %±^, (5.8) 

A) 

D,ix; A) D {x; t(A + (i - 1)5)) (5.9) 

^e{x; A) 
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Compare with the similar expressions for the Xi Laguerre & Jacobi polynomials (2.10)-(2.15) 
& (3.13)-(3.18) in [H], and for the Xe Wilson & Askey- Wilson polynomials (3.7)-(3.9) in 

m. 



Since det^^(A) = ni=cr y Bt{x;X) ^ for the parameter range under consideration, 
the operators ^^(A) and Ai{\y have no zero modes. By using the two formulas fl4.6l) - fl4.7p . 
we can show that 

(A) = k,{X) [mx + iS + S) + /,,o(A)6,,o(A)) , (5.10) 
•H^)(A) = ke{X){'Hi{X) + /,,o(A)6,,o(A)), (5.11) 



where ^^(A) is 



Therefore the original system with the shifted parameters ('H(A + £5 + 6)) and the deformed 
system ('H^(A)) are exactly isospectral. Note that the maximal value of x for ^/(A + id + S) 
is N — i{= x'^^^). Based on the results fl5.10l) - fl5.11l) . we have 



(+)^x;A) = 0„(x;A + ^(5 + 5), ^['J {x; X) = (j)e,n{x; X) , (5.13) 



4t (^) = l^e{X){Sn{X + iS + S) + /,,o(A)6,,o(A)) = /«,(A) (f,,„(A) + /,,o(A)6,,o(A)). (5.14) 

The correspondence of the pair of eigenfunctions 0^^^(x) with their own normalisation spec- 
ified in the preceding sections are related by 



■^(A)0j+)(x; A) 
x; A) = A).,(A) ^7 (5.15) 



i,(A)t&Ha;; A) 



Ul;X)sfXX) Vf^i{X)be^n{X) 



Let us introduce the operators ^^(A) = {^g.j.^y{X)) and Be{X) = {Bf-^^y^X)) {x,y = 0,1, . . . 
a^max) defined by 



/•,(A) = a/6(1;A)s,(A)^,(x; X)-' o ^^i^ o 0o(a;; X + iS + 6), (5.17) 

a/kHA) 

i3,(A) = , ^ 0o(x; A + M + ^)-i o 4= o M^; A). (5.18) 

6(l;A)s,(A) V'^HA) 
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Their explicit forms are: 

^KA) = - , ^ ^} A + - 1)6) - (x; A + (£ - 1)5)6"^) . 

^e{x;X)ip{x;X+{i-l)X + 5)^ ' 

(5.20) 

Compare with the similar expressions for the Xn Wilson & Askey- Wilson polynomials (3.20)- 
(3.21) in [15]. The operators ^^(A) and i3f(A) act as the forward and backward shift operators 
connecting the original polynomials P„ and the exceptional polynomials P(,^n- 

T,{\)Pn{x- A + M + 5) = /£,„(A)P,,„(x; A), (5.21) 
E,[\)P,^n{x; A) = 6,,„(A)P„(a;; A + £(5 + I). (5.22) 

The former relation ( ]5.2ip with the explicit form of J-'^(A) ( 15.19^ provides the explicit ex- 
pression (14.281) of the exceptional orthogonal polynomials. In terms of -^^(A) and i3^(A), the 
relations ( ]5.10p -( l5.1ip become 

Bi{\)T(,{X) = -HiX + i6+~d) + /,,o(A)6,,o(A), (5.23) 
J-,(A)i3,(A) = -HeiX) + /,,o(A)6,,o(A). (5.24) 

The other simple consequences of these relations are 

4?(A) = k,{X)knWk,nW, Sn{\ + id) = A„(A)6,,„(A) - /,,o( A)6,,o (A) . (5.25) 

The i'^ inner product for and (pe^m can be calculated in the following way: 



■; A),0£,m(-;A)) 



1 lUl; A)s,(A) ^^^^^^ _ . A), i,(A)0™( ■ ; A + M + 5)) 



= T^A ^^%7Xr^ (A(A) V...( ■ ; A), 0.( . ; A + £5 + ^)) 
= -H^6(l; A)s,(A) (0„( ■ ; A + M + 5), 0™( ■ ; A + M + ^)) 

/£,m(A) 

= ^,(1; A) "^""^ , ^f'"^^^ s,(A), (5.26) 

'rf„(A + M + (5)2/.,n(A) 

where we have used (15.131) . (I5.16P and (12.81) . This gives a proof of (I4.40p . 
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It is interesting to note that the operator ^^(A) intertwines those of the original and 
deformed systems ^(A) and ^^(A): 

AiiX + S)A{X + iS + S)= Ae{\)AeiX), (5.27) 
^(A)^(A + M + Sy = Ae{XyAe{X + S). (5.28) 

In terms of the definitions of the forward shift operators J^(A) fl2.18p . J^i{X) (14.331) . ^^(A) 
flSTTj) . and B{X) (127791) . Bi{X) fOD . the above relations are rewritten as: 

se{X + S)Pe{X + 6)T{X + £6 + 6) = se{X)Te{X)fi{X), (5.29) 

se{X):Fe{X)B{X + £6 + 6) = se{X + 6)Be{X)Pe{X + 6), (5.30) 

where Si{X) is 

MA)'=^«:.(A)x| + (5.31) 

These relations can be proven by explicit calculation with the help of the two formulas of 
the deforming polynomial ^^(x; A) fl4.6l) - fl4.7l) . 

By applying Ai{X + 6) and ^^(A) to (1232|) and (l213ll (with replacement A ^ A + M + ^) 
respectively, together with the use of (15.271) . (15.281) and (15.161) . we obtain 



Ae{X)(j)i^n{x] X) 



lk e{X + 6) sejX) Ul;X) 1 kn-i{X + 6) 

keiX) seiX + 6)^,{l;X + 6)B{0;X + £6 + ~6) /,,„,(A) 
X /„(A + M + 6)^i,n-i{x; X + 6) 

^ /„(A + M)0,,„_i(x;A + (5), (5.32) 



'5,(0; A) 



A / Ke{X) se{X + 6)^e{l;X + 6) ^ , ,x , x\ /^."W 

Ae{X)^^e,n-i{x;X + 6) = \ . , ttt ^———B{0;X + £6 + 6)^ - 

" Ke{X + 6) se{X) U^;X) /,,„„i(A + <5) 

X 6„_i(A + £6 + 6)(j)i^n{x; X) 



= y/BfXO; X) 6„_i(A + £6)(j)e,n{x; X + 6). (5.33) 

In the calculation use is made of the exphcit forms of ke{X), ^^(A), Bi{x] A), /,^„(A), /n(A) 
and fen (A) in the second equalities. This provides a proof of (I4.30p - p.32p without recourse 
to the shape invariance of the deformed system. Likewise the above intertwining relations of 
the forward-backward shift operators (I5.29p - (l5.30p give a proof of (I4.36p - (l4.37p . respectively, 
again without recourse to the shape invariance. 
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Since the g-Racah polynomial P^{x] A) f l3.10p is related to the Askey- Wilson polynomial 

p„(cos x; a, 6, c, as |[5j 

Pf{x;X)= . " . p^{l{d-2q- + d-h-n-,ad--2,bd-"^,cd~-^,d-2\q), (5.34) 
(a,6,c;g)„ 

many formulas for the g-Racah case in sections S] and are obtained essentially from those 
for the Askey- Wilson case [13] by the following replacement: 

gixAw ^ ^A^w ^ (5.35) 

6 Other Polynomials: dual (g')-Hahn, little Q'-Jacobi 

In §|3]-§|5] we have derived the exceptional Racah and g-Racah Hamiltonian systems by 
deforming those of the Racah and g-Racah in parallel in terms of a degree i polynomial 
with twisted parameters. It is well known that the Racah polynomials can be obtained from 
the g-Racah polynomials by taking the standard g — )■ 1 limit with an appropriate overall 
rescaling. The same limiting procedure could be applied to derive the exceptional Racah 
polynomials from the exceptional g-Racah polynomials. 

Likewise various orthogonal polynomials of a discrete variable can be obtained from the q- 
Racah polynomials by many different limiting procedures with/without the g — )■ 1 limit. Here 
we present two such examples: the dual (g)-Hahn and the little g-Jacobi polynomials and the 
corresponding exceptional polynomials. The former is a finite dimensional example and the 
latter is infinite dimensional. It should be stressed, however, that there is no guarantee that 
the limiting procedure among the undeformed polynomials could be lifted to produce the 
corresponding exceptional polynomials. For example, the Hermite polynomials are known to 
be obtained from the Jacobi or the Laguerre polynomials by a certain limit procedure. But 
that does not produce exceptional Hermite polynomials from the known exceptional Jacobi 
or Laguerre polynomials. 

6.1 Dual (g')-Hahn 

In this subsection we present the ordinary and the exceptional dual Hahn (dH) and the dual 
g-Hahn (dgH) polynomials. Like as (g)-Racah cases, these are finite dimensional: Xmax = 
'^■max = N and x^^^^ = n^^^ = N — £. The dual g-Hahn case is obtained from the g-Racah 
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case by the following limit: 

q^"^ = {q~^ ,a,t,abq~^), qR dgH. 



(6.1) 



The dual Hahn case is obtained from the dual g-Hahn case by taking g — )■ 1 hmit with an 
appropriate overall rescaling. 

6.1.1 Original systems 

The Hamiltonian systems thus obtained belong to the e = 1 case of pLj and they are listed 
as follows: 



X={a,b,N) : dH au a u 

q^ = {a,b,q^) :dgH ' ^ = M : dR,dq}i, n 

a > 0, 6 > : dH 

0<a<l, 0<6<1, : dgH ' 

{x + a){x + a + b - 1){N - x) 



1 : dH 
1 : dgH 



B{x;X) 




dH 
dgH 



{2x- 1 + a + b){2x + a + b) 
^qX-N _ l)(l _ ag^)(l - a6g^-i; 

(1 -a6g2^-i)(l -a6g2^) 
x{x + b - l){x + a + b + N - 1) 
{2x-2 + a + b){2x -1 + a + b) 

(1 - g")(l - a6g"+^"i)(l - bq^-'] 



(l-a6g2'^-2)(l -a6g2^-i; 



dH 

dgH 



: dH 
1 :dgH 



r]{x; X) 

302 



x{x + a + b — 1) 
(g-^' - 1)(1 - abq'^-^] 
-n, X + a + b — 1, —X 

a, -N 
'g~", abq^~^, q~ 

a, g-^ 

Rr,iri{x;X);a-l,b-l,N) : dH 

i?„(l + abq'^ + 7]{x] X) ; ag~\ 6g"\ N\q) : dgH ' 
[a)x {2x + a + b — l){a + b) 



q;q 



: dH 
:dgH ' 

: dH 
: dgH 



x\{N-xy. {b)^{x + a + b-l)N+i 

{q;q)N {a,abq~^ ■,q)x 1 - abq^''-'^ 

{q ', <l)x {q ; q)N-x {abq^, b ; q)^ a^' 1 - abq'^ 



(a)„ {b)N-n (&)iV 
— X 



dniXf 



n\{N-n)\ {b)N " {a + b)N 

{q;q)N {a;q)n{b;q)N-n ^ ib]q)Na 
{q;q)n{q]q)N-n {b]q)Na'' {ab^q)^ 



N 



dH 
dgH 
dH 
dgH 



(6.2) 
(6.3) 

(6.4) 

(6.5) 
(6.6) 



(6.7) 



(6.8) 



(6.9) 
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^{x- X) 



2x + a + h 

a + b 
q^^ — abq^ 

1 — ah 

6.1.2 Deformed systems 



dH 



: dgH 



/„(A)=^„(A), bniX) = l :dH,dgH. 



(6.10) 



We restrict the parameter range of (16. 3 p as follows: 

a > 0, 6 > 1 : dH ffi in 

< a < 1, < 6 < g, : dgH ■ ^ ^ ^ 

The data for the Hamiltonian systems of the exceptional dual (g)-Hahn polynomials are as 

follows: 

= Pe{x;i{X+{i-l)d)), t(A) = (Ai + A2 + A3-l,l-A3,l-A2) : dHdgH 
■-i, a + b + x + £-2, -X 



3-^2 



3^2 



a + b + N -l,b-l 
q;q 



q ^, abq^^^ ^, g ^ 



vf{x-X) 



abq^-\ bq-^ 

(x - N){x + a) 

a + b-l 
.,.(l-g--^)(l-ag-) 



: dH 
: dgH 



dH 
dgH 



(x; A) 



v^{x-X) 



<(x;A) 



1 — abq ^ 

X + a + b -1 
; : dH 

a + b-l 

_^1 - abq""-^ ' 

g : dgH 

1 — abq ^ 

{x + a + b + N -l){x + b-l) 

a + b-l 

^_.^_i^i_^(l^a6g-'+^-i)(l-6g-i: 



X 



1 — abq ^ 
: dH 



: dH 
: dgH 



a + b-l 

abq^^- T~Zi '■ ^Q'H 



1 — abq 
S = (0,-1,0) : dH,dgH, 

-b-N + n + l 
_rigi-^(l-6g^-"-i 



ke{X) = 



: dH 
: dgH 



, 6£,„,(A) = 1 :dH,dgH, 



1 : dH 

6g^-^-i :dgH 



(1-6) 



a+b+N -I 



: dH 



a+b+i-l 



(6.12) 
(6.13) 
(6.14) 

(6.15) 

(6.16) 

(6.17) 
(6.18) 

(6.19) 
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s^(A) = ki{X) : dH,dgH. (6.20) 
Note that /^„(A), s^(A) < and bi^nW > 0. All the formulas in §[3]-§[5]are satisfied. 

6.2 Little g-Jacobi 

In this subsection we present the ordinary and the exceptional little g-Jacobi {IqJ) poly- 
nomials. They are infinite dimensional: Xmax = '"-max = oo and x^ax = ''^max = The 
Hamiltonian system of the little g-Jacobi polynomials is obtained from that of the g-Racah 
polynomials by the following limit: 

q^""" = {q~^,aq''^h-\bq,t-'), qR ^ algH IgJ, (6.21) 

where algH stands for the alternative g-Hahn system (with A = {aq,bq, N)) in §5.3.1 of [1]. 

6.2.1 Original system 

The data of the shape invariant Hamiltonian system whose eigenf unctions are described by 
the little g-Jacobi polynomials are as follows [1]: 

q^ = {a,b), (5 = (1,1), K = q-\ < a < q-\ < b < q-\ (6.22) 

5(s;A) = a(g-^-6g), D{x] X) = q'"" - 1, (6.23) 

^„(A) = (g-"-l)(l-a6g"+i), r]ix; X) = 1 - q^ , (6.24) 
P.(x;A) = P„(r/(x;A);A) = 30i(^ ",a6gp,g ^ | ^ . ^- y j 

= (-a) g 2 ^ 'j7 — ^2011 g;g ) 

{bq;q)n \ aq J 

= (_a)-V^"^"+'^7^^Pn(l-r](x;A);a,6|g), (6.25) 
[bq ; g)n 



^>o(x;A)2 = ^^^(ag)^ (6.26) 

{q;q)x 



2 (bq, abq ; q)nd'^q'^ 1 — abq'^'"'^^ {ciq]q)c 



(g, aq ; g)„ 1 - a6g [abq^ ; g)oo 

</.(x;A) = g^ /„(A)=£„(A), 6„(A) = 1. (6.28) 

6.2.2 Deformed system 

The data for the exceptional little g-Jacobi polynomials are as follows: 
6(x; X) = Uv{x;X+{^-m-A) 
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= Pt(x-, t(A + (C - 1)«)), t(A) S (-Ai - 2, A2) 



<(a;; A) = -ag^+\ (x; A) = 1 - (6.30) 

t;f (x; A) = -g^ (x; A) = 1 - (6.31) 

^ = (1,-1), (6.32) 

/An(A) = g-"(l-ag"+i)^-^, 6,,„(A) = 1 - (6.33) 

k,{\) = {aq'+')-\ s,(A) = ^^, (6.34) 

s,{X) = k,{\){l-bq'). (6.35) 



Note that /^^„(A), 6£,„(A), Sf(A) > 0. All the formulas in §[3]-§[5]are satisfied. 

7 Summary and Comments 

The Racah and the g-Racah polynomials are the most generic members of the orthogonal 
polynomials of a discrete variable satisfying second order difference equations. By deforming 
the discrete quantum mechanical systems governing these polynomials in terms of degree 
£ eigenpolynomials, the exceptional Racah and g-Racah polynomials are obtained as the 
main part of eigenfunctions of the deformed systems, which are shape invariant and exactly 
solvable. By certain limiting procedures, the exceptional dual (g)-Hahn polynomials and the 
exceptional little g-Jacobi polynomials are derived. The deformation process goes parallel 
with that for the exceptional Wilson and Askey- Wilson polynomials. Some of the charac- 
teristics of the quantum mechanics with real shifts are the cause of complications which led 
to the delayed discovery. The method of deriving the exceptional polynomials is new to the 
theory of orthogonal polynomials. As for the parameter ranges in which the orthogonality 
weight functions are positive, we have made a quite conservative arguments. It is quite 
possible that for a fixed i the valid parameter range could be enlarged than those given in 
the text. On the other hand, the difference equations for the original and the exceptional 
orthogonal polynomials, (I2:22|) - fl2:25| ). (iOejl - dCTj) and ^^-^^ are purely algebraic 
and they hold for any parameter values. 

With the understanding of all the generic exceptional orthogonal polynomials as solu- 
tions of exactly solvable quantum mechanical systems, the next challenge would be the 
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construction of the exceptionals of various reduced cases, for example, the Morse potential, 
the Meixner-Pollaczek and the Krawtchouk cases, etc. Finding multivariable generalisation 
is truly interesting but its feasibility is as yet unclear. 
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